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Abstract 

We employ a self-consistent simulation approach based on quantum theory to investigate the physical 
properties of a pair of ferromagnetic and antiferromagnetic nanotubes. It was observed that under the 
given conditions, no matter the external magnetic field was absent or applied along the easy longitudinal 
axis, the spins always ordered in that direction due to the special geometric shape of the tubes and the 
magnetic uniaxial anisotropy, so that the two sorts of nanosystems exhibit typical ferromagnetic and 
antiferromagnetic properties, which may find applications in modern technology, in strong contrast to 
the phenomena observed previously in nanoparticle [T], where an external magnetic field, applied parallel 
to the antiferromganetically coupled spins, is able to turn the spins off their original direction to form 
symmetric pattern around the field direction. This peculiar feature enable us to build up one-dimensional 
chain models for the two sorts of nanosystems. Considering their fast computational speed and simplicity, 
these theoretical models were then utilized to investigate the finite size effects of the nanosystems, and 
perform further analysis for long tubes. Especially, our results obtained with the theoretical models and the 
numerical approach are exactly identical, verifying the correctness and applicability of the computational 
methodology. 
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Since the discovery of carbon nanotubes [2], the field of this kind of nanosystems, especially 
those exhibiting magnetic properties, have been attracting considerable attention due to their 
potential applications in modern technology as well as their significance in science [11S1 [7]. 
For instance, high attention has been paid to the magnetic nanotubes based on the transition 
metals, such as CoPt, CoPd, FePt and FePd alloys owing to their potential applications in high- 
density magnetic recording media 0 m [H HU [H]. Now, one can find various new nanotubes 
based electronic and spintronic devices m, multi-functional systems involving nanotubes for 
biotechnological applications M- 

So far, many techniques for fabricating magnetic nanotubes have been developed [muni Hum 
Hiiin]- Hematite a-Fe 2 03 nanotubes can be produced for example by using a template method, 
which provides nanotubes with diameters between 100 nm and 200 nm, and a wall thickness 
less than 10 nm [16]; attachment of maghemite 7 -Fe 203 nanoparticles to the surface of carbon 
nanotubes via a modified sol-gel technique has also produced nano hybrid structures behaving as 
nanotubes m- 

Theoretically, electronic structures of ferromagnetic single-wall nanotubes were investigated 
by means of spin density functional theory EU; continuous models for ferromagnetism [22] . ap¬ 
proaches based on Green functions [23] , and simulation methods like Monte Carlo [24] , have been 
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also employed to investigate the exotic magnetic properties and spin configurations of magnetic 
nanotubes. 

In the present work, we first investigate the magnetic properties of a pair of ferromagnetic and 
antiferromagnetic nanocylinders, in which Heisenberg and uniaxial anisotropy along the central 
longitudinal z-axis co-exist, by means of a self-consistent approach based on quantum mean field 
theory. We find that the special geometric shape of the tube surfaces and the uniaxial anisotropy 
force the spins order in the easy z direction, so that the two sorts of nanosystems exhibit typical 
ferromagnetic and antiferromagnetic properties, which may find applications in modern technology, 
in strong contrast to the phenomena observed in antiferromagntic nanoparticle [T], where an 
external magnetic field, parallel to the oppositely aligning spins, is able to turn the spins off their 
original direction to form symmetric pattern around the field direction. This peculiar feature 
of the nanosystems enables us to build up one-dimensional ferromagnetic and antiferromagnetic 
chain models for theoretical analysis. It turns out that the magnetizations, hysteresis curves, etc., 
generated by means of the numerical approach and the theoretical models are exactly identical, 
verifying our computational methodology once again. 

1 Computational Model 

Now, we consider a nano-cylinder, which is formed by rolling up a rectangular monolayer lattice 
consisting of S =1 spins sitting on the square lattice sites. The length of the nanotube is £ = iV^a 
along the longitudinal central axis, where a is the lattice parameter, and there are JVji spins 
on every circle around the axis, i.e., the circumference of such a circle is C = JVfia. Thus the 
coordinators of a spin on the nanotube wall are determined by 
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^ 2TTn' 

Vn 

= — sm 
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( 1 ) 


where 1 < n < Nji, and 1 < m < N^. 

To do Monte Carlo simulations for Ising-like nanocylinders, Masrour et. al. only considered 
Heisenberg exchange interaction among the nearest neighboring spins with uniform strength, the 
uniaxial anisotropy along the longitudinal central axis, and the external magnetic field exerted in 
that direction [23]. Following Masrour’s work, the Hamiltonian for our single wall nanotubes thus 
can be expressed as 


H E E - ffSMB E 


B 


* j/* 


( 2 ) 


where gs=2 is the Land factor, Jij and Ka represent the strengths of the exchange interaction 
among the neighboring spins and the uniaxial anisotropy assumed to be in the z-direction as well. 
In the above Hamiltonian, the spins are quantum operators instead of the classical vectors. Since 
5 = 1, the matrices of the three spin components in Cartesian coordinate system are given by 


^/ 0 \/ 20 \ / 0 \/ 20 \ / 100 \ 

& = - C2 0 ^/2 , Sy = j-A -V2 0 -V2 , Sz = 0 0 0 

^ \ 0 V2 0 J \ 0 V2 0 J \ 0 0 -1 J 


(3) 


respectively. 

Since only the exchange interactions between the nearest spins are considered, a spin on the 
tube wall interacts only with its two nearest neighbors in the same line parallel to z axis with a 
strength of J \, two nearest neighbors on the same circle around the central axis with a strength of 
Ji. When both Ji and Ji are positive, the nanotube is ferromagnetic, the spins tend to align in 
the same direction; while both Ji and Ji are negative, the nanotube is antiferromagnetic, every 
pair of neighboring spins attempt to order in the opposite directions. 
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This Hamiltonian given by Eq.Q was first used to study the magnetic properties of the nan¬ 
otubes numerically by means of the SCA approach [TJ 1251126) . In every simulation step, it was 
diagonalized so as to calculate the magnetic moment of the considered spin. All of our simula¬ 
tions were started here from a random magnetic configuration and a temperature well above the 
magnetic transition, then carried out stepwise down to very low temperatures with a temperature 
step AT < 0. At any temperature, if the difference |(^') — {Si)\ between two successive iterations 
for every spin is less than a very small given value tq, convergency is considered to be reached 

2 Results Obtained Numerically with the SCA Approach 





Figure 1. Simulated magnetic structures projected onto (a) j/ = 0, and (b) ^ = 15o cross sections, respectively; 
calculated {Sz) (c) as function of temperature, and (d) as the function of « at different temperatures, for the 
ferromagnetic nanotube by means of the SCA approach in the absence of external magnetic field. Here Nl = 30. 

Firstly, we performed simulations for a ferromagnetic nanotube with = 30 and Nn = 20 by 
means of the SCA approach. To do this, Heisenberg exchange and uniaxial anisotropy strengths 
were assigned to J7i = 1 K, ^2 = 0.8 K and = 0.1 K, respectively. Fig.l(a,b) displays the 
magnetic structures projected onto the y = 0 and z = 15a cross sections. As shown there, all spin 
order ferromagnetically antiparallel to the z-axis below Curie temperature Tq « 2.44 K. That 
is, in the magnetic phase, only {Sz) is nonzero, both {Sx) and (Sy) all vanish. As displayed in 
Fig.1(c), while temperature arises, (Sz) decreases gradually from the maximum value at very low 
temperatures to zero until Tc- Due to the finite length of the nanotube, we can find from Fig.1(d) 
that the spins near the two ends of the nanotube naturally have weaker magnitudes than those in 
the middle part at all recorded temperatures. 
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Figure 2. Simulated magnetic structures projected onto (a) y = Q, and (b) z = 15a cross sections; calculated (Sz) 
(c) as the function of temperature, and (d) as the function of 2 ; at different temperatures, for the antiferromagnetic 
nanotube by means of the SCA approach in the absence of external magnetic field. Here Nl = 30. 


For comparison with the results just obtained, we carried out simulations for an antiferromag¬ 
netic nanotube of the same size scale with the SCA approach by only changing the signs of the 
Heisenberg exchange constants, that is, now = -1 K and J\ = -0.8 K. Fig.2(a,b) depict the 
magnetic configurations projected onto the y = 0 and 2 = 15a cross sections. There we hnd that 
all spin order antiferromagnetically in the 2 -direction below Nel temperature T/v ~ 2.44 K. Once 
again, in the magnetic phase, only (Sz) is nonzero, both (Sx) and (Sy) all vanish. As temperature 
arises, the two opposite (S'z)’s, depicted in Fig.2(c), attenuate gradually from the saturated value 
S' =1 at very low temperatures to zero at Tat. As Nz approaches to 1 or from the middle 
region, |(Sz)| decreases as seen in Fig.2(d). However, these curves exhibit excellent symmetry 
both horizontally and vertically. 

The total free energy F, total energy E, magnetic entropy Sm and specihc heat Cm of these 
canonical systems can be evaluated with 

F =-kBT log Zn, E = —-^logZff, 

Sm = + kslogZN, CM=T{^^)g , (4) 

successively, where /3 = l/^ksT) and Zn is the partition function of the whole system. Figure 
3(a,b) display the F, E, and Cm curves obtained with the SCA approach for the antiferromagnetic 
nanotube. It is very interesting to find that for the counterparting ferromagnetic tube, i.e., as l7i 
= 1 K and J 2 = 0.8 K, the calculated F, E and Cm curves are almost exactly identical. The 
sudden change of E and the sharp peak in the Cm curve near Tjv are the signs of phase transition. 
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Figure 3. (a) Total energy and total free energy, and (b) specific heat, per mole of spins calculated with the 
SCA approach for the antiferromagnetic nanotube in the absence of external magnetic field. Here Nl = 30. 




Figure 4. Hysteresis curves of the (a) ferromagnetic, and (b) antiferromagnetic nanotubes calculated with the 
SCA approach at different temperatures. Here Nl = 30. 


To learn how these two sorts of nanosystems react to the external magnetic field, we then 
studied their hysteresis processes with the SCA approach. Fig.4(a,b) show the longitudinal hys¬ 
teresis curves calculated at four different temperatures below Tm with the SCA approach for the 
ferromagnetic and antiferromagnetic nanotubes, respectively. In the ferromagnetic case, hysteresis 
loops are formed at all four temperatures. Especially at T = 0.4 K, the loop is almost a rectangle. 
If the external magnetic field is applied antiparallel to the original magnetization, all spins are 
suddenly rotated for 180° at B = Bc= 1.25 T for instance. In the antiferromagnetic case, no 
hysteresis loop can be formed as shown in Fig.4(b). And in the both cases, (Sx) and (Sy) always 
vanish if the magnetic field is exerted exactly in the z-direction. In our previous SCA simulations 
for an antiferromagnetic nanoparticle, it was discovered there that a strong external magnetic 
field, applied in the direction along which the antiferromagnetically coupled spins ordered inside 
the core region, could rotate them off their original direction, forming symmetric pattern around 
the field direction [T]. The geometrical shape of the nanotubes has indeed endowed them with 
very typical ferromagnetic and antiferromagnetic characters. 
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3 One dimensional ferromagnetic and antiferromagnetic chain models 


As just described, when a weak uniaxial anisotropy along the longitudinal axis is present, all spins 
order spontaneously in the easy z-direction. Even when an external magnetic field is applied along 
the central longitudinal axis, nonzero (Sx) and {Sy) could not be observed. These peculiar features 
allow us to build up one dimensional ferromagnetic and antiferromagnetic chain models. In the 
ferromagnetic case, there are Nl equally spaced spins on the line along the longitudinal axis. The 
m-th spin Szim) interacts with its left neighbor Szim — 1) and its right neighbor Sz(jn + 1) with 
strength J7i, other two neighbors Sz(jn) on the same circle around the central axis with strength 
Ji- Since the spins only order in the z-direction, we then have, for the Hamiltonian given in 
Eq.(I5|), three pure eigenvectors = | + 1), |■^/) 2 ) = |0), and l^a) = | — 1) for the m-th spin, and 
the corresponding eigenvalues are expressed as 


£1 = - \J\ {{Sz{m - 1)) + {Sz{m + 1))) + 2 J 2 {Sz{m)) + gsHsB] - Ka , 

£2 = 0 , (5) 

£3 = [Ji {{Szim - 1)) + {Sz{m + 1))) + 2J2{Sz{m)) + gsgsB] - Ka , 


respectively. This chain has a finite length N^. However, if we let 0 < m < + 1, (S'z(O)) 

and {Sz{Nl + 1)) assigned to 0 at the beginning, thus all spins can be treated equally. Further 
assuming ^ = Ji {{Sz{m — 1)) + {Sz{m + 1))) + 2J2{Sz{m)) + gsgsB, we finally get a formula 


2 exp {KA/ksT) sinh j^/ksT) 

1 + 2 exp (ATA/fcsT) cosh ’ 


( 6 ) 


by making use of quantum theory. 

An antiferromagnet can be considered to be composed of two oppositely oriented A and B 
subsystems. Thus, in the simplified model, an A type spin interacts only with the nearest H-type 
spins, and vice versa. More specifically, an m-th A spin Sf{m) interacts with a B-type neighbor 
Sf{m — 1) on its left, another B-type neighbor Sf{m -b 1) on its right, both with strength Ji, 
and other two B type neighbors S^{m) belonging to the same circle around the central axis with 
strength ^72, and vice versa. Thus, for the same sake just described, we have once again three 
pure eigenvectors = 1 + 1)) = |0)) = I “ 1) 1®^ Hi® m-th A or B spin, 

respectively, and the corresponding eigenvalues are given by 

ef= - [Ji {{Sf’^im - 1)) + +(m + 1))) + 2J 2 (Bf+(m)) + gs/rsB] - Ka , 

= 0 , (7) 

[^i((Bf+(m-l)) + (Bf+(m + l)))+2J2(^i'’^(m))+gs/rBB] 

respectively. Assuming gA,B = Ji {{Sf’^{m - 1)) -b {Sf’^{m + 1))) -b 2 J 2 {Sf’^{m)) + gsgsB, 
and making use of quantum theory, one can easily deduce a pair of coupled formulas 


2 exp {KA/ksT) sinh {ijB/kBT) 

1 -b 2exp{KA/kBT) cosh (Tys/fesT) 
2 exp {KA/ksT) sinh {gA/ksT) 

1 + 2exp {KA/kBT) cosh {gA/kBT) 


for the two sorts of oppositely oriented spins. 


4 Results obtained with the one dimensional magnetic chain models 

Using the one-dimensional chain models, we have successfully reproduced the results given in 
Fig.(1,2,4), no matter if the external magnetic field was absent or applied in the z direction. 
There is not need to present those results here once again. Obviously, computational speed can 
be greatly accelerated by using these theoretical models because of their simplicity in comparison 
with the SCA approach. 
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T(K) 



Figure 5. Calculated (Sz) (a) as the functions of temperature for the ferromagnetic nanotubes with different 
lengths, (b) as the functions of at different temperatures for a ferromagnetic tube of length Nl = 10, by means 
of the one-dimensional ferromagnetic chain model in the absence of external magnetic field. 


With the two theoretical models, it is much easy to focus our studies on the finite size effects. 
If the periodic condition is applied, we then have an infinitely long nanotube. Fig.5(a) depicts 
the spontaneous magnetization curves calculated with the one-dimensional chain model for three 
ferromagnetic nanotubes of finite and infinite lengths, respectively. We find here that a short tube 
has a low Curie temperature Tq, however the length does not strongly affect the Tq value. To 
fully understand the size effects, we plot the magnetization curves for the tube with a length of 
= 10 at three different temperatures in Fig.5(b). Near the two ends, the spins there have weaker 
magnitudes than those in the middle region; and at higher temperatures, these effects become 
much stronger. Fig.5(b) looks very similar to Fig.1(d). But now, because of the short length, the 
size effects, for example at T = 2 K-i spread deeply into the very middle part of the nanotube. 
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Figure 6. Calculated (Fz) (a) as the functions of temperature for the antiferromagnetic nanotubes with different 
lengths, (b) as the functions of W at different temperatures for an antiferromagnetic tube with Nl = 10, by 
means of the one-dimensional antiferromagnetic chain model in the absence of external magnetic field. 


For comparison, the magnetization curves for the three antiferromagnetic tubes, as displayed 
in Fig.6, were calculated with the one-dimensional chain model by only changing the signs of J\ 
and Ji used in Fig(5) . Now we find the similar size effects as just described for the ferromagnetic 
nanotubes. And more interestingly, if we only draw the upper parts of Fig.6(a,b), the curves will 
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Figure 7. Calculated hysteresis curves of the (a) ferromagnetic, and (b) antiferromagnetic nanotubes at different 
temperatures by means of the one-dimensional ferromagnetic and antiferromagnetic chain models. Here Nl = 


200 . 


coincide exactly with their counterparts depicted in Fig.5(a,b) correspondingly. Therefore, we can 
say the two sorts of systems are symmetric, and our calculated results are self-consistent. 

The relatively faster computational speed of the theoretical models makes it easy for us to 
do simulations for long nanotubues. Fig.7(a,b) show the longitudinal hysteresis curves of the 
ferromagnetic and antiferromagnetic nanotubes of a length =200, obtained by means of the 
theoretical models using the parameters given herein. The curves shown here for the long nan¬ 
otubes actually differ slightly from those shown in Fig.4(a,b), even though their lengths differ 
four times. Once again, the nanosystems exhibit typical ferromagnetic and antiferromagnetic 
characters. 
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Figure 8. Calculated longitudinal susceptibility curves for (a) a ferromagnetic, and (b) an antiferromagnetic 
nanotubes with different lengths, by means of the one-dimensional ferromagnetic and antiferromagnetic chain 
models, respectively. 


Finally, we calculated the susceptibilities of the two sorts of nanosystems. For the purpose, a 
very weak external magnetic field is usually considered to be exerted in a special direction, and the 
susceptibility can be approximately estimated as the ratio of the induced/enhanced magnetization 
to the applied magnetic field. So, we defined the longitudinal susceptibility as Xc ~ {Sz)/B, and a 
field of 0.01 Tesla was assumed to be applied along the z-axis to calculate the quantities with the 


























5 Conclusions and Discussion 


9 


one-dimensional magnetic chain models. The unit of the Xc is 1/(Tesla Per spin). For simplicity, 
it is omitted in the figures. 

Fig.8(a) displays the calculated results for three ferromagnetic nanotubes of different lengths. 
At very low temperatures, (Sz) is close to the saturated value 5=1, thus Xc ss {Sz)/B fv 
100/(Tesla Per spin). As temperature arises, Xc attenuates. The applied magnetic field modify 
each curve’s shape around Tc considerably so that they attenuate continuously. Due to the size 
effects, we can observe here three distinct curves. 

In the antiferromagnetic case, as T arises below Tjv, the applied magnetic field gradually rotate 
those magnetic moments antiparallel to it to its direction. So in Fig.8(b), we observe that Xc 
increases until Tat, where a peak appears. But above Tjq, the disordering effects of temperature 
becomes stronger and stronger gradually, so we find all Xc’s fade with increasing temperature. 
Again, we can see three curves near Tat’s, which is the evidence of size effects. 

5 Conclusions and Discussion 

We have also performed simulations by exchanging the magnitudes of Ji and 172, that is, by 
assigning Ji = ±0.8 K and 172 = ± 1 K, respectively, but we obtained the same results as those 
shown in Fig.(1-4). This is easy to understand. For instance, when l7i = 0.8 K and 172 = 1 K, a 
spin still interacts with two spins with a strength of 1 K, other two spins with a strength of 0.8 
K, as in the case when l7i = 1 K and Ji = 0.8 K. So the calculated results with the two sets of 
parameters are naturally same. 

In summary, we firstly used the SCA approach to simulate the magnetic structures for the 
two sorts of magnetic nanotubes, and calculated their spontaneous magnetization, longitudinal 
hysteresis curves and thermodynamic quantities. We found that no matter the external magnetic 
field was absent, or applied along the central longitudinal axis, the spins always aligned in that 
direction, seeming all bundled to the tubes’ surface because of their peculiar geometric structure, 
so that the tubes exhibit typical ferromagnetic or antiferromagnetic features, which may find 
applications in industry. In contrast, an external magnetic field, applied in the direction along 
which the antiferromagnetically coupled spins align inside a nanoparticle, is able to turn the 
spins off the line to form symmetric pattern around the field direction [1]. Above findings in 
numerical simulations enable us to build up one-dimensional ferromagnetic and antiferromagnetic 
chain models, which could reproduce all the results obtained in SCA simulations in principle, 
and were then employed to investigate the size effects, and investigate the magnetic properties 
of long tubes. Because of the simplicity of the theoretical models, computational speed can be 
considerably improved. Especially, the exact agreements between our numerical calculations and 
theoretical analysis justify the correctness of our quantum simulation approach once again. 
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